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TIME  AND  LOAD  INDEPENDENT  PROPERTIES 
OF  BITUMINOUS  MIXTURES 

SYNOPSIS 

This  paper  presents  a  new  approach  to  characterize  the  time -dependent 
behavior  of  asphaltic  concrete  by  a  transfer  function  which  is  derived 
from  the  frequency  response  of  the  material,  using  sinusoidal  loading  tests. 
The  ratio  of  the  Laplace  transform  of  the  output  to  the  Laplace  transform 
of  the  input  to  any  system  is  defined  as  the  transfer  function  of  the 
system  between  the  input  and  the  output.   The  transfer  function  is  a 
function  of  the  complex  variable  s  ■  jw  ,  where  to  is  the  frequency  of  the 
£nput  and  j  is  the  imaginary  unit. 

In  this  study  asphaltic  concrete  was  treated  as  a  damped,  linear  system 
in  which  the  applied  force  f(t)  was  the  input  and  the  resulting  displacement 
x(t)  was  the  output.   Both  the  input  and  output  were  functions  of  time.   By 
varying  the  frequency  of  the  input  the  output- input  magnitude  ratio  was 
obtained  for  each  frequency  and  its  plot  against  log  frequency  yielded  the 
frequency  spectrum.  Also,  the  time  lag  between  the  output  and  the  input  was 
noted  from  which  the  phase  angle  was  calculated. 

A  simple  geometric  technique  to  obtain  the  transfer  function  directly 
from  the  frequency  spectrum  was  used  in  this  investigation.   This  technique 
was  applied  to  the  results  of  the  test  specimens  of  asphaltic  concrete  cut 
from  laboratory  compacted  beams  of  two  different  compositions  tested  at  three 
temperatures.   It  is  shown  that  the  transfer  function  for  bituminous  concrete 
thus  derived  from  a  dynamic  test  can  be  used  to  predict  the  displacement  of 
the  test  specimen  subjected  to  a  static  load  by  treating  the  latter  as  a  step 
function  of  time  and  through  the  use  of  Laplace  inverse  transform.   The 
excellent  agreement  between  the  calculated  and  measured  values  of  the  dis- 
placements indicates  that  the  transfer  function  represents  a  material  property 
which  is  independent  of  load  input. 


2. 
The  concept  of  transfer  function  makes  it  possible  to  represent  the  time- 
dependent  behavior  of  asphaltic  concrete  by  a  fourth  order  linear  differential 
equation  with  constant  coefficients  without  assuming  any  spr ing-dashpot  model. 
The  coefficients  can  be  computed  from  the  roots  of  the  denominator  and 
numerator  of  the  transfer  function. 

INTRODUCTION 

Theoretical  developments  in  the  stress-displacement  behavior  of  asphaltic 
concrete  in  a  pavement  structure  are  very  restricted  in  their  applicability 
to  represent  actual  field  conditions  due  to  the  complexity  of  asphaltic 
concrete  as  a  structural  material  and  the  complicated  interaction  with  the 
other  structural  layers.   A  rigorous  theoretical  solution  to  a  typical 
boundary  value  problem  in  pavement  mechanics  must  satisfy  the  equations  of 
equilibrium  (or  motion),  compatibility  equations  (or  some  other  form  of 
conservation  equation),  boundary  conditions,  and  initial  conditions.   The 
development  of  such  a  solution  requires  the  use  of  some  form  of  stress- 
displacement-time  relation  for  the  material,  and  it  is  this  aspect  which 
specifies  the  solution  for  a  particular  medium.   At  present,  there  are  no 
satisfactory  constitutive  relations  available  for  pavement  materials,  both 
asphaltic  concrete  and  soils,  and  this  is  perhaps  the  greatest  impediment  to 
realistic  theoretical  solutions  of  the  response  of  pavement  systems  to  applied 
loads,  dynamic  or  static. 

Having  the  ultimate  objective  to  design  such  layered  pavement  structures, 
the  theory  of  linear  viscoelasticity  has  been  widely  used  for  the  evaluation 
of  pavement  components  as  well  as  for  understanding  their  response  to  varied 
loading  conditions.   Layered  systems  were  not,  however,  analyzed  using  time- 
dependent  (viscoelastic)  material  properties  until  after  the  development  of 
the  correspondence  principle  for  isotropic  media  in  1955  by  Lee  (l)*  and  shortly 


*  Figures  in  parentheses  refer  to  references  in  the  Bibliography. 


thereafter  the  extension  by  Biot  (2)  to  include  anisotropic  media.  Stress 
and  strain  for  linear  viscoelastic  materials  can  be  related  by  either 
differential  or  integral  linear  operators.  The  differential  operator  form 
of  the  stress-strain  law  is  most  commonly  used  and  as  Biot  (3;  has  pointed 
out  may  be  visualized  as  a  combination  of  springs  and  dashpots.  More  general 
methods  of  viscoelastic  stress  analysis  have  been  described  by  Lee  {k)  and 
Lee  and  Rogers  (5). 

Pister  and  Monismith  (6)  have  applied  the  basic  differential  equation 
for  the  three  element  model  to  the  solution  of  a  viscoelastic  beam  on  an 
elastic  foundation,,   Fister  (7)  has  considered  the  solution  of  a  viscoelastic 
plate  en  a  viscoelastic  foundation.   Other  literature  (8,,  9?  10)  also  exist 
in  limited  quantity  dealing  with  viscoelastic  slabs  or  foundations  in 
various  ways.  Harr  (11)  used  a  two  element  model  to  show  the  influence  of 
vehicle  speed  on  pavement  deflection.  The  usefulness  of  these  methods  is 
restricted  to  the  extent  of  the  realistic  representation  of  the  viscoelastic 
characteristics  of  the  materials  of  the  pavement  structure  in  the  numerical 
analysis. 

The  study  of  viscoelastic  characteristics  of  bituminous  mixtures  has 
been  the  object  of  many  research  workers,  and  here  again  spring-dashpot 
models  have  been  generally  made  use  of.   By  using  a  Burger's  model.  Wood 
and  Goetz  (12)  found  that  the  behavior  of  a  sheet  asphalt  mixture  obeyed 
the  laws  of  linear  viscoelasticity  for  limited  stresses  and  small  deformations, 
Secor  and  Monismith  (13)  observed  that  the  four  element  model,  suggested  by 
Kuhn  and  Rigden  (lU)  was  capable  of  representing  the  displacement  of  an 
asphaltic  concrete  only  to  a  limited  extent. 

Dynamic  tests,  such  as  those  employing  sinusoidal  loading,  have  been 
used  to  study  the  viscoelastic  response  of  asphaltic  concrete,  by  some  investi- 
gators (15,  16,  IT) .  They  have  attempted  to  characterize  the  viscoelastic 
response  by  using  the.  phase  angle  and  the  complex  elastic  modulus,  obtained 


from  the  input-output  traces  of  the  dynamic  tests.   However,  these  parameters 
are  not  easily  adaptable  for  numerical  analysis.   Besides,  they  are  functions 
of  the  input  load,  which  seriously  limits  the  ability  of  these  parameters  to 
truly  represent  the  load  response  characteristics  of  the  material. 

This  motivated  the  authors  to  search  for  a  parameter  or  parameters  or  a 
function  which  can  represent  the  viscoelastic  response  of  the  material  under 
any  load  and  which,  at  the  same  time,  is  independent  of  the  type  and  magni- 
tude of  the  input  load.   Such  a  function  should  necessarily  be  time-dependent 
and  should  be  used  as  the  material  characteristic  for  both  static  and 
dynamic  loads.   A  study  of  the  literature  reveals  that  transfer  functions  are 
frequently  used  to  characterize  a  dynamic  system  in  electrical  and  mechanical 
engineering  problems  ( 18  through  23)-   While  serving  as  a  link  between  the 
time -dependent  input  and  output  of  the  system,  the  transfer  function  completely 
describes  its  dynamic  response.   Hence,  it  was  attempted  to  investigate  the 
feasibility  of  applying  the  concept  of  transfer  functions  to  describe  the 
response  of  bituminous  mixtures  to  dynamic  and  static  loads.   This  paper 
presents  a  method  of  obtaining  the  transfer  function  for  a  given  bituminous 
concrete  from  laboratory  tests  and  studies  the  extent  to  which  this  function 
can  represent  the  viscoelastic  characteristics  of  the  material.   Effects 
of  mix  type,  temperature  and  anisotropy  on  the  transfer  function  were  studied 
with  a  view  to  providing  insight  into  this  function  which,  it  is  hoped,  may 
facilitate  the  development  of  more  refined  methods  for  the  structural  design 
of  pavements. 

It  is  assumed  in  this  investigation  that  asphaltic  concrete  behaves  as 
a  linear  system.  Although  it  has  been  recognized  that  the  behavior  of  this 
material  is  nonlinear  particularly  at  higher  stress  levels  and  at  higher 
temperatures,  the  investigations  of  Busching,  Goetz  and  Harr  (2^),  and 
others  (15,  25)  have  established  the  validity  of  this  assumption  for  small 
displacements  normally  encountered  in  bituminous  mixtures. 


CONCEPT  OF  TRANSFER  FUNCTION 
The  ratio  of  an  operational  output  (Laplace  transform  of  the  output)  of 
a  dynamical  system  to  the  operational  input  (Laplace  transform  of  the  input) 
is  called  the  transfer  function  between  the  operational  input  and  its 
corresponding  output.  The  transfer  function  is  shown  schematically  in  Fig.  1 
It  is  by  definition  a  function  of  the  complex  variable  s. 


Operational 
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Fig.  1.   Block  diagram  representation 
of  the  transfer  function. 


If  a  rigid  body  is  subjected  to  a  dynamical  force  f(t),  then  the  force 
and  the  resulting  displacement  x(t)  can  be  considered  as  the  input  and  out- 
put, respectively  for  the  dynamic  system.  The  transfer  function  between  the 
operational  force  and  the  operational  displacement  is  given  by 

G(s)  «  iLli  ...       (!) 

£(s) 

where     G(s)  -  the  transfer  function 

x(s)  -  Laplace  transform  of  the  output,  w  x(t)> 
f(s)  =  Laplace  transform  of  the  input,  i  |_f(t)}. 

Solving  for  x(s), 

x(s)  -  G(s)f(s)  ...       (2) 

The  inverse  transform  of  x(s)  is  x(t)  and  is 

x(t)  »  r1{G(s)f(s)j  ...       (3) 


6. 
Eq.  3  shows  that  once  the  transfer  function  G(s)  is  known  for  any  system, 

the  displacement  x(t)  can  be  evaluated  for  that  system  for  another  given 

input  force  f(t). 

In  general,  the  transfer  function  G(s)  is  a  ratio  of  two  polynomials 

in  s,  such  as 

N(s)     (s  "  ax)(s  '  a2)  •••  (s  '  an) 
G(s)=#f  =   (s  -  bl)(s  -  b2)  ...  (s  -  bj 

where  a   's  and  b.  's  are  the  roots  or  zeros  of  the  numerator  N(s)  and 
denominator  D(s)  of  G(s)  respectively.   b.  's  are  also  called  the  poles  of 
G(s),  since  the  transfer  function  becomes  infinite  when  evaluated  at  a  zero 
of  its  denominator.   That  is, 


G(b.)  =  "  (5) 

Substituting  for  G(s)  in  Eq .  2,  x(s)  can  be  rewritten  as 


f(s)(s  -  a1)(s  -  a2)  ...  (s  -  an) 
X(S)  =  (s  -  bl)(s  -  b2)  ...  (s  -  bj  (6) 

Eq.  6  can  be  expanded  in  m  partial  fractions  for  the  m  roots  of  D(s)  as 


which  is  an  identity  in  s  where  C,  ,  C^,  ....  C   are  constants  independent  of 

l7      2  m 

s.      Using   the   Laplace    inverse   transform,    the    final    solution   for   x(t)    is 
obtained    from  Eq .    7  as 

b   t  b   t  \t 

x(t)=C1e1        +C2e2       +..-+Cme  (8) 

From  the  above  discussion  it  is  obvious  that  once  the  transfer  function 
for  a  system  is  known,  the  displacement  x(t)  of  the  system  for  a  given  force 


input  f(t)  can  be  evaluated  theoretically. 

DETERMINATION  OF  THE  TRANSFER  FUNCTION 
The  experimental  basis  to  determine  the  transfer  function  for  a  test 
specimen  is  the  frequency  spectrum,,  which,,  in  turn,  may  be  obtained  from 
sinusoidal  input-output  data  of  the  test  specimen  (23) »  The  plot  of  the 
amplification  (that  is,,  the  ratio  of  the  output  magnitude  to  the  input 
magnitude)  against  the  test  frequency  is  the  frequency  spectrum^  and  for 
convenience  the  magnitude  of  the  amplification  is  plotted  in  decibels**. 
The  frequency  spectrum  is  approximated  with  a  series  of  connected  straight 
lines  so  as  to  form  asymptotes  to  the  curve,,  see  Fig.  2.  The  intersections 
of  the  asymptotes  determine  the  corner  frequencies.  The  modified  transfer 
function  is  then  obtained  (26)   as  a  function  of  the  slopes  of  the  asymptotes 
and  their  corner  frequencies  as 

1        2  1        3  2  n  n-1 

G^jw)  «  A(j<*  u^)6  (p+    i£)  6   (joj+u)3)  6  ...  (ju+wn) 

The  modified  transfer  function  is  the  vector  amplification  or  simply  the 

ratio  of  the  magnitude  of  the  output  to  that  of  the  input  and  is  a  function 

of  the  input  frequency.,   By  substituting  s  -  j  w  in  Eq.  9?  the  transfer 

function  for  the  system  is  obtained 


G(s)  -  A  (s  +w1;  t 

The  constant  A  can  easily  be  determined  from  any  one  experimental  frequency 
in  the  frequency  spectrum  by  calculating  the  absolute  value  of  the  transfer 
function  from  Eq,  9   at  that  frequency  and  equating  it  to  the  corresponding 
amplification  value  in  the  frequency  spectrum. 


*  Rdb  =  20  logio  R 
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9* 
In  approximating  the  frequency  spectrum,  the  slopes  of  the  asymptotes 

are  generally  taken  to  be  six  decibels  per  octave*  for  greater  ease  of 

determination^,  although  straight  lines  of  any  slope  may  be  used.  The  closer 

the  approximation  to  the  actual  frequency  spectrum  curve^  the  better  will 

be  the  approximation  of  the  transfer  function  being  sought „ 

EXPERIMENTAL  INVESTIGATION 

The  experimental  phase  of  this  study  had  as  an  objective  the  develop- 
ment and  use  of  accurate  techniques  for  obtaining  all  of  the  data  required 
in  the  determination  of  the  transfer  function  from  a  frequency  spectrum. 
The  core  cf  the  technique  lies  in  devising  a  system  which  can  apply  a 
sinusoidal  force  input  of  a  desired  magnitude  and  frequency  and  which  can 
simultaneously  measure  the  displacement  output „ 

An  MIS  electronic  function  generator  coupled  with  a  loading  frame  which 
was  fitted  with  an  electronically  controlled  hydraulic  actuator  was  used 
for  this  p"Tpose„   See  Figs*  3  and  k*     The  loading  frame  was  installed  in 
a  constant  temperature  room  which  facilitated  carrying  out  the  tests  at 
the  desired  temperatures.  A  two-channel  Erush  recorder  which  was  attached 
to  the  function  generator  recorded  the  input  and  output  of  the  test  specimen 
simultaneously . 

Dynamic  Tests 
The  function  generator  was  used  to  apply  sinusoidal  force  inputs  to  the 
asphaltic  concrete  specimen  by  suitable  manipulation  of  the  controls.   The 
magnitudes  of  the  force  varied  from  2.5  to  15  lbs.  and  the  frequencies  varied 
from  0o00l6  to  1.6  cycles  per  second. 


*  An  octave  is  defined  as  the  frequency  range  for  which  the  ratio  of  the 

upper  %->und  frequency  to  the.  lower  bound  frequency  is  2  to  1;  that  is, 

a)  <w  2u>.   is  an  octave. 
1      1 
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11. 

Since  in  a  sine  test  the  specimen  is  alternatively  in  compression  and 
in  tension,,  the  method  cf  mounting  the  specimen  in  the  test  assembly  was 
of  paramount  importance.   Two  aluminum  plates,  2  in,  x  2  in.  x  3/8  in.  in 
size j,  were  glued  to  the  specimen  top  and  bottom  respectively,  using  commer- 
cially available  Eastman  910  Adhesive.  The  plates,  in  turn,  were  connected 
to  the  top  and  bottom  platens  by  screws  and  nuts.   The  top  platen  was 
rigidly  fixed  to  the  moving  ram  of  the  hydraulic  actuator  while  the  bottom 
platan  was  fixed  to  the  loading  deck.   Sec  Fig.  5°   Remotely  controlled  and 
commanded  by  the  electronic  console  of  the  function  generator  j,  the  hydraulic 
ram  moves  up  and  down  at  the  pre-set  sine  force  and  frequency.  Thus  the 
specimen  was  subjected  to  the  desired  sinusoidal  force  input  which  was 
recorded  by  the  left  channel  of  the  Brush  recorder.  The  corresponding 
displacement  of  the  specimen  was  sensed  by  the  L7DT  in  the  actuator  and  was 
simultaneously  recorded  as  the  output  in  the  right  channel  of  the  recorder. 
Fig.  6  shows  the  two  traces  thus  obtained  in  a  typical  dynamic  test. 

Static  Tests 
The  same  test  assembly  which  was  used  for  the  dynamic  tests  was  used 
for  the  static  tests.  When  the  function  generator  is  set  for  a  combination 
of  &   ramp  function  and  a  high  frequency^,  the  actuator  applies  a  static 
compressive  force  of  desired  magnitude.  As  before,  the  left  and  right 
channels  recorded  the  input  and  output  respectively.  The  two  traces  of  a 
typical  static  compression  test  are  illustrated  in  Fig.  7° 

Materials  and  Preparations  of  Specimens 
Two  different  gradings.  one  with  a  maximum  size  of  the  No.  k   sieve  and 
the  other  the  3/8  inch  sieve  were  used  in  this  investigation.  The  bituminous 
mixes  prepared  on  the  basis  of  these  gradings  were  designated  as  Mix-1  and 
Mix-2_,  respectively.  A  60-70  penetration  asphalt  cement  was  used.  The  mix 
was  compacted  into  beams  of  size  2  in.  x  2  1/2  in.  x  12  in.  from  which  test 
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Fig.  5*   Specimen  in  position  for  testing 


13 


Fig.  6.   Typical  traces  of  a  sinusoidal  loading  test 
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Fig.  7.  Typical  traces  of  a  static  compression  test 


15- 
specimens  of  dimensions  1  in.  x  1  in.  x  2  in.  were  cut  along  the  x,  y  and  z 

axes  of  the  beam,  x  being  the  longest  axis  and  g  th€  ?;rtical  axis,   The 

specimens  were  designated  as  X-l,  Y-l,  Z-l,  X-2,  Y-2  and  Z-2,  the  letter 

denoting  the  axis  along  which  it  was  cut  and  the  numeral  the  mix,   The  details 

of  obtaining  the  specimens  are  shown  in  Reference  26. 

DISCUSSION  OF  RESULTS 

Frequency  Response  of  Asphalt ic  Concrete 
In  most  of  the  dynamic  tests  a  magnitude  of  10  lbs.  was  employed  for 
the  sinusoidal  load  input  for  frequencies  covering  three  decades,  varying 
from  0.01  radians  per  second  to  10.0  radians  per  second.   These  were  observed 
to  be,  respectively,  the  slowest  and  the  fastest  frequencies  to  which  the 
test  specimens  were  responsive.  The  input  magnitude  and  the  output  dis- 
placement were  recorded  in  pounds  and  inches,  respectively.  Also,  both  of 
them  were  recorded  in  millivolts. 

The  ratio  of  the  output  to  the  input  expressed  in  decibels  (also  referred 
to  as  amplification  or  gain)  plotted  against  leg  frequency  gives  the  desired 
frequency  spectrum.  The  phase  lag  or  simply  the  phase  angle  for  each  fre- 
quency was  also  computed  from 

0  *=(*)  t  (11) 

where  0  -   the  phase  angle  in  radians 

w  =  the  test  frequency  in  rad/sec,  and 
t  -  time  lag  between  input  and  output  in  seconds. 
The  phase  angle  is  also  plotted  against  log  frequency.   The  two  plots  together 
represent  the  frequency  response  of  the  test  specimen  completely. 

Table  1  shows  typical  test  results  for  a  dynamic  test  and  Fig.  8  gives 
the  frequency  response  for  the  same  test.  It  is  seen  that  the  displacement 
of  the  specimen  is  largest  at  the  slowest  frequency  and  is  least  at  the 
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TABLE  1 
Results  of  A  Typical  Sinusoidal  Test 


Specimen: 


z-1 


Test  Temperature;   80  F 


CO 

rad/sec 


1 

rnv 


? 
mv 


R   =  — 
Sl 


Nib 


t 

sec 


rad     0     deg, 


0.025 

0.050 

0.100 

0.250 
0.500 

1.000 

2.500 
5.000 

10.000 


1000 


22.0 
19.5 

17.0 
11.0 

9.0 

8.0 
6.5 
5.0 

k.o 


0.0220 
0.0195 

0.0170 
0.0110 
0.0090 

0.0080 
0.0065 
0.0050 

0.00^0 


•33.1520 

.■$k.k2k0 

-35-3920 
-39° 1720 
40.9160 

41.9380 

43.7^20 
46.0200 

47.9580 


15-30 
8.32 

it.  15 

1.50 

0,68 


0.382 
0.kl6 

0A15 

0.375 

0.3^0 


0.311  0.311 

0.110  0.271 

o.o46  0.230 

0.021  0.210 


21.8 
23.9 

23.8 
21.5 

19.5 

17.8 

15.5 
13.2 

12.0 


R      - 


LEGEND 
input  frequency,  rad/sec 
input  magnitude^  millivolts 
output  magnitude^  millivolts 


a2/a1 


db 


■  R  expressed  in  decibels 


20  loglcR 


t  =  time,  lag  between  input  and  output,  sec, 
0     =  phase,  lag  between  input  and  output 
=  00  t  radians  . 
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fastest  frequency.   Thus  the  gain  continuously  falls  with  increase  in 

frequency  which  indicates  the  overdamped  nature  of  the  test  specimen. 

The  phase  angle  increases  with  increase  in  frequency  in  the  first 

decade  of  test  and  then  decreases  so  that  a  bell-shaped  curve  results  when 

phase  angle  is  plotted  against  frequency.   The  peak  has,  in  general,  been 

observed  to  occur  in  the  region  of  0.05  to  0.1  rad/sec.   These  results 

substantiate  the  observations  of  Pagen  and  others  who  have  studied  the 

dynamic  response  of  bituminous  concrete  in  sinusoidal  testing  (15 ,    16,  17)- 

Form  of  the  Transfer  Function 
Table  2  shows  a  set  of  typical  transfer  functions  obtained  for  the 
specimens  tested  at  80F,  from  which  it  is  seen  that  the  transfer  function 
derived  for  each  individual  case  is  of  the  form 

A(s  +  ai)(s  +  ag)(s  +  aj 

G(S)  =   (s  +b1)(s  +hz){s   +  b^)(s  +bk)  (12) 

where  A  is  a  constant  and  a.  's  and  b.  's  are  the  corner  frequencies  as 

11 

previously  defined  and  which  are  obtained  in  the  process  of  approximating 
the  frequency  spectrum  with  asymptotes.   The  significant  feature  of  this 
form  of  the  transfer  function  is  that  the  factors  in  both  the  numerator 
and  denominator  are  to  first  power.   It  may  be  recalled  that  this  is  due  to 
the  geometric  procedure  followed  in  approximating  the  frequency  spectrum 
with  asymptotes  of  slopes  of  six  decibels  per  octave.   It  was  also  mentioned 
that  any  slope  can  be  used,  in  which  case  it  is  obvious  from  Eq.  10  that 
the  transfer  function  may  contain  fraction  powered  or  higher  powered  terms 
in  either  the  numerator  or  the  denominator  or  both.   This  form  is  not 
desirable  since  the  resulting  transfer  function  will  be  too  difficult  or 
even  impossible  to  handle  through  Laplace  transforms. 

From  the  procedure  followed  in  approximating  the  frequency  spectrum 
with  asymptotes  it  is  obvious  that  the  number  of  factors  in  the  numerator 


TABLE 


Transfer  Functions  at 
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Specimen 
z-1 

x-1 

y-1 

z-2 


x-2 


y-2 


G(s] 
G(s] 
G(s] 
G(s] 
G(s] 
G(s; 


Transfer  Function 

(0.05)(s  +  0.025)(s  + 

0.2)(s  +  2) 

(s  +  O.G2)(s  +  0.1)(s 

+  l)(s  +  10) 

(0.0lH5)(s  +  0„028)(s 

+  0.17)(s  +  1.8) 

(s  +  0.02)(s  +  0.l)(s 

+  1)(  s  +  10) 

(0.04)(s  +  0.025)(s  + 

0.l6)(s  +  2.3) 

(s  +  0.02)(s  +  0.l)(s 

+  l)(s  +  10) 

(O.OU)(s  +  Oo025)(s  + 

o.i6)(s  +  2.3) 

(s  +  0.02)(s  +  0.l)(s 

+  l)(s  +  10) 

(0.0Ul5)(s  +  0.028)(s 

+  0.17)(s  +  1.8) 

(s  +  0.02)(s  +  0.l)(s 

+  l)(s  +  10) 

(0.0Ul5)(s  +  O.C28)(s 

+  0.17)(s  +  1.8) 

(s  +  0.02)(s  +  0.l)(s 

+  l)(s  +  10) 
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and  denominator  of  the  transfer  function  in  Eq.  12  depends  on  the  number  of 

asymptotes  used  to  approximate  the  frequency  spectrum.   It  was  found  from 
the  analysis  of  the  experimental  results  that  the  transfer  function  obtained 
using  eight  asymptotes  was  accurate  enough  for  practical  purposes. 

Effect  of  Mix  Type  on  Transfer  Function 
The  two  different  bituminous  mixes  studied  in  this  investigation  were 
quite  dissimilar  in  their  aggregate  gradings  and  in  their  binder  contents, 
namely,  S.Ofo  and  3-5$  DY  weight  of  aggregate  for  Mix-1  and  Mix-2,  respect- 
ively.  However,  it  is  apparent  from  Fig.  9  that  there  was  no  appreciable 
difference  in  the  transfer  functions  of  the  specimens,  irrespective  of  their 
composition  or  orientation  when  tested  at  any  one  temperature.   Even  though 
the  grading  and  asphalt  content  were  different  for  the  two  mixes,  their 
density-void  ratio  characteristics  were  similar,  namely,  142.6  pcf  and 
6.%   for  Mix-1  and  142.0  pcf  and  Q.%   for  Mix-2. 

It  can  be  anticipated  that  density  will  have  a  marked  influence  on  the 
response  to  loading  of  a  given  mix  and  hence  on  its  transfer  function. 
However,  more  work  is  necessary  to  study  the  effect  of  this  and  other  mix 
variables . 

Effect  of  Specimen  Orientation 
It  was  pointed  out  under  "Experimental  Investigation",  that  for  each 
mix  three  specimens  were  cut  along  three  mutual  perpendicular  directions, 
designated  x,  y  and  z.   It  is  seen  from  Fig.  9  that  there  is  not  much 
difference  in  the  transfer  functions  of  the  differently  oriented  specimens 
of  either  mix  at  any  one  temperature.  This  suggests  that  effect  of 
anisotropy,  if  any,  is  not  reflected  in  the  transfer  function. 

Effect  of  Temperature 
The  experiments  in  the  present  investigation  were  conducted  at  70, 
80  and  90F .   By  examining  the  absolute  value  curves  in  Fig.  9,    it  is 
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readily  seen  that  temperature  plays  an  important  role  in  the  transfer 

function  of  any  specimen.   Asphalt  cement  being  thermoplastic,  the 
asphaltic  concrete  specimen  will  yield  increased  displacements  for 
increases  in  temperature  at  a  given  stress  level.   This  results  in  a 
decrease  in  viscosity  of  the  binder  in  the  asphaltic  concrete  as  temper- 
ature increases.   Thus  the  output-input  ratio  for  a  constant  input  at  any 
frequency  will  increase  with  decrease  in  viscosity,  or  increase  in  temper- 
ature.  This  is  precisely  what  is  shown  by  the  curves. 

Examination  of  the  transfer  functions  revealed  that,  for  a  given 
specimen,  the  various  factors  in  the  numerator  and  denominator  do  not 
seem  to  vary  much  with  change  in  temperature,  but  that  the  coefficient 
A  varies  appreciably.   It  would  be  of  interest  to  study  this  change  in  A 
with  reference  to  the  viscosity  of  asphalt  cement  or  the  asphalt  cement- 
filler  matrix. 

Effect  of  Specimen  Size  on  Transfer  Function 
Frequency-response  tests  were  conducted  on  three  different-sized 
specimens  with  constant  height-width  ratio  using  Mix-1  at  90F.   They  were 
cut  along  the  same  direction  in  the  compacted  beam.   The  large,  medium 
and  small  sizes  were  1  l/k"   x  1  1/V  x  2  1/2",  1"  x  1"  x  2"  and  3/V  x 
3/V1  x  1  1/2",  respectively.   Analysis  of  their  test  results  showed  that 
the  frequency  spectrums  for  the  three  specimens  are  nearly  the  same  and 
that  they  can  be  approximated  by  one  single  transfer  function,  which  shows 
that  the  transfer  function  for  the  mix  is  independent  of  the  specimen  size. 

In  the  light  of  this  observation,  it  appears  that  the  transfer  function 
can  be  used  to  represent  the  dynamic  characteristics  of  a  viscoelastic 
material  in  much  the  same  way  as  the  elastic  modulus  represents  the  stress- 
strain  characteristics  of  an  elastic  material.   In  other  words,  for  a  given 
viscoelastic  material  such  as  a  bituminous  concrete  there  is  only  one 
transfer  function  for  a  given  temperature. 


Prediction  of  Displacements  for  Sinusoidal  Loads 
It  was  seen  under  "Concept  of  Transfer  Function"  that  once  the  trans- 
fer function  of  a  system  is  known,,  it  can  be  used  to  predict  the  output 
of  the  system  for  any  other  given  input  which  is  a  function  of  time,  After 
the  transfer  functions  for  the  asphaltic  concrete  test  specimens  under 
investigation  were  determined,  attempts  were  made  to  predict  the  output 
displacement  of  the  specimens  when  subjected  to  a  sinusoidal  input  of  known 
magnitude,  other  than  the  one  used  to  determine  the  transfer  function. 
Each  specimen  was  tested  under  three  sinusoidal  inputs  at  each  of  the  three 
temperatures  used  previously. 

From  Eq.  3*  the  solution  for  the  output  displacement  of  a  dynamic 
system  with  a  transfer  function  G(s)  for  a  sinusoidal  load  inpat  of  f(t)  « 

f  sin  at  is  given  by 

o  J 

x(t)  ■  f  |  G(s)  I  sinfot  +  0)  (13) 

where  f  ■  magnitude  of  the  input 
0)   =  frequency  of  the  input 
0  =  phase  angle  between  output  and  input. 
In  Eq.  13 }   x(t)  is  a  maximum  when  sin(a)t  +  0)  ~  1„  so  that 

Experimentally,  it  is  convenient  to  measure  x(t)    at  the  peak  of  the 
sinusoidal  displacement  output.  Thus,  the  calculated  displacement  from 
Eq.  lit-  can  be  compared  with  the  measured  displacement  for  a  given  magnitude 
and  frequency. 

It  was  mentioned  elsewhere  that  in  the  experimental  set-up  of  this 
investigation,  the  input  was  recorded  in  pounds  and  millivolts  and  the 
output  was  recorded  in  inches  and  millivolts.  Thus  the  output- input  ratio 
can  be  either  dimentionless  or  in  units  of  inch  per  pound.  The  frequency 
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spectrum  and  the  transfer  function  analyses  were  obtained  as  dimensionless 

quantities.   In  applying  units  to  Eq.  lk ,   an  experimental  constant  K  needs 

to  be  introduced  and  it  becomes 

x(t)    =  K  f  |G(s)|  (15) 

v  'max      o    '  ' 

From  the  sensitivity  control  values  of  the  recorder,  K  was  calculated  and 
found  to  be 

-k 
K  =  25  x  10   inch/pound  (l6) 

The  calculated  and  measured  values  for  the  maximum  displacement  for 
all  the  specimens  tested  at  80F  are  shown  in  Table  3-   In  all  cases  the 
magnitude  of  the  sinusoidal  input  was  kept  the  same  and  only  the  frequency 
was  varied.   Comparison  of  the  calculated  and  the  measured  displacement 
values  clearly  indicates  the  close  agreement  between  these  values  in  all 
cases.   Besides  proving  the  efficacy  of  the  transfer  function  as  a  displace- 
ment predicting  tool  for  viscoelastic  materials,  this  also  shows  that  the 
bituminous  concrete  behaves  as  a  linear  system  at  the  levels  of  stress 
considered  in  the  tests. 

Prediction  of  Displacements  for  Static  Loads 
The  transfer  function  can  be  used  to  predict  the  time-dependent  dis- 
placements of  a  system  subjected  to  a  step  function  input,  and,  mathemati- 
cally, a  static  load  can  be  represented  by  a  step  function.   For  the 
transfer  function  given  by  Eq .  12,  the  solution  for  the  displacement  under 
a  static  load  can  be  shown  as 


b  t      b  t      b  t      b.t 
x(t)  =  re    +  Ce    +CeJ  +  C,  e  H  +C  (17) 


where  C   's  are  constants  independent  of  t  for  the  system  and  b.  's  are 
i  l 

the  roots  of  the  denominator  of  the  transfer  function  which  are  observed 
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TABLE  3 
Calculated  and  Measured  Displacements  at  80  F 


Specimen 

Sine 

Load  Input 

Constant  Load 

Input 

f 
o 

lbs 

to 

rad/sec 

x(t) 

f 
0 

lbs 

x! 

O 

calculated 
inch  (xlQ  ) 

measured 
incn(xI0  ; 

calculated 
inch(xlO   ) 

measured 
inch(xlO*4) 

z-1 

5 
5 
5 

0.05 
0.5C 
5-00 

2.37 
1.20 
O.58 

2.25 
1.25 
0 .  5C 

5 

:: 
15 

3-13 
6.26 
9-39 

2.5C 

5-00 
10.00 

x-1 

5 

5 
5 

0.05 
0.50 
5-00 

1.59 
0.8? 
0.50 

1.25 
0.75 
0.50 

5 
10 

15 

2.2^ 

6.69 

2. CO 
3.75 
5-50 

y-1 

15 
15 
15 

0.05 
0.50 
5-00 

5-32 
3ol8 
1.43 

5.00 
3-00 
1.75 

5 
10 

15 

2.30 

k.6o 

6.90 

2.25 
4,00 

8.75 

z-2 

15 
15 
15 

0.05 
0.50 
5-00 

5-32 
3-18 
1.43 

5.00 
3-00 
1.50 

5 

10 
15 

2.30 
4.60 
6.90 

2.25 

4.50 
6.50 

x-2 

15 

15 
15 

0.05 
0.50 
5-00 

4.76 

2.62 
1.50 

5.00 
2.50 
1.75 

5 
10 

15 

2.23 
4.46 

6.69 

1.25 

4.00 

9.50 

y-2 

15 

15 
15 

0.05 
0.50 
5.00 

I+.76 
2.62 
1.50 

^.75 
2.50 
1.50 

5 
10 

15 

2.23 
k.k6 

6.69 

2.25 
4.00 
7.00 
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to  be  negative  in  the  case  of  the  experimentally  derived  transfer  functions 

Two  boundary  conditions  are  applicable  to  Eq .  17:   l)  when  t  =  0,  and 

2)  when  t  =  ». 


When  t  =  0, 


that  is 


x(t)  =  0  (16) 


C5  =  Cl  +C2  +C3  +  Ck  (19) 

When  t  =  % 

x(t)  =  C5  (20) 

When  a  constant  load  is  applied  to  a  viscoelastic  material,  the  dis- 
placement approaches  a  constant  value  after  a  certain  time  depending  upon 
the  nature  of  the  material.   In  the  laboratory,  for  the  bituminous  concrete 
specimens  tested,  this  time  was  observed  to  be  in  the  order  of  a  few  minutes. 
Hence,  the  second  boundary  condition  discussed  above  can  be  applied  to  this 
steady  displacement  value  in  the  test,  that  is,  Eq .  20  can  be  used  to  cal- 
culate the  displacement.   The  limitations  of  using  Eq .  20  are  recognized 
in  that  in  a  test  t  does  not  reach  infinity  and  hence  the  measured  value  will 
be  slightly  lower  than  the  calculated  value. 

In  order  to  get  the  calculated  displacement  in  proper  units,  Eq .  20 
should  be  multiplied  by  the  experimental  constant  K,  described  in  the  previous 
section  so  that 

x(t)  =KC5  (21) 

where  K  =  25  x  10   inch/pound.   C   has  the  units  of  pound. 

The  calculated  and  measured  values  of  the  displacement  at  80F  for  the 
specimens,  each  under  three  different  constant  loads,  are  shown  in  Table  3- 
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Comparison  of  the  calculated  and  the  measured  displacement  values  clearly 

indicates  the  close  agreement  between  the  two  in  almost  all  cases.  This 
observation  is  very  significant  and  it  brings  to  focus  the  following  points: 

1.  The  mathematical  theory  of  transfer  functions  is  applicable  to 
viscoelastic  materials  in  general,  and  to  bituminous  concrete  in  particular. 

2.  The  techniques  developed  to  derive  the  transfer  function  in  this 
investigation  are  sound. 

3.  The  transfer  function  serves  as  a  connecting  link  between  the 
responses  of  the  material  tested  under  dynamic  and  static  loads. 

Differential  Equation  from  Transfer  Function 
Once  the  transfer  function  is  known  for  a  given  system,  the  different- 
ial equation  for  the  system  behavior  can  be  written  using  the  Laplace  inverse 
transform.  For  bituminous  concrete,  the  transfer  function  was  seen  to  be 
of  form 


6(s) 


A(s  +  a1)(s  +  ag)(s  +  a^) 
(s  +  bJTfs  +  bg)(s  +  b^)(s  +  bk) 


(12) 


From  definition 


From  Eqs.  1  and  12, 

f(s) 
Rewriting, 


G(s)  -  ^^ 
f(s) 


A(s  +  a  )(s  +  a  )(s  +  a  ) 


(s  +  bx)(s  +  b  )(s  +  b  )(s  +  b^) 


(1) 


(22) 


£(s)(s  +  bx)(s  +  b2)(s  +  b3)(s  +  b^)  =  A  f(s)(s  +  a^s  +  a2)(s  +  a3)   (23) 

Multiplying  and  rearranging  the  terms, 

x(s)[s  +  B1s3  +  B2s2  +  B  s  +  B J  »  A  f(s)[js3  +  B  s2  +  BgS  +  B  ]  (24) 
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where   B  =  b  +  b  +  b  +  b, 

B2  =  bxb2  +  b2b3  +  b^  +  b^  +  b^  +  b3bl 
B3  -  blb2b3  +  b2b3b^  +  b3b1+bi  +  b^blb2 

\  "   blb2b3\ 

B5  =  ax  +  a2  4-  a3 

B6  =  3la2  +  a2a3  +  a^ 
B7  =  ala2a3 

Rewriting  Eq .  2k, 

k-  ■?-         2- 

s  x(s)  +  B  sJx(s)  +  B  s  x(s)  +  B  sx(s)  +  B,x(s) 

=  Afs3f(s)  +  B  s2f(s)  +  BgSf(s)  +  B  f(s)J 

Applying   the   Laplace    inverse    transform, 

"^  +  Bi   Am  +  B2  Afcl  +  .  *fcl  +  v(t) 


(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 


(32) 


dt 


dt 


=  A 


r^  &j|l  *fcl  +      (t) 

dt  at 


(33) 


Eq-  33  is  a  fourth  order  linear  differential  equation  with  constant  coeffi- 
cients which  describes  the  dynamics  of  the  bituminous  concrete.   The  input 
(force)  and  the  output  (displacement)  are  the  two  variables  in  the  equation 
which  are  functiorS  )f  time.   The  constants  B.  's  are  easily  determined  from 
the  roots  of  the  denominator,  b.  's,  and  from  those  of  the  numerator,  a.  's, 
of  the  transfer  function,  using  Eqs.  25  through  31- 

It  is  of  significance  to  note  here  that  these  roots  of  the  differential 
equation  are  obtained  as  the  corner  frequencies  directly  from  the  frequency 
spectrum  for  the  material  under  test,  and  no  assumption  whatsoever  has  been 
made  in  deriving  the  equation. 
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CONCLUSIONS 

Based  on  the  results  and  within  the  limitations  of  this  investigation, 
the  following  conclusions  are  enumerated? 

1.  The  viscoelastic  or  time -dependent  characteristics  of  an  asphaltic 
concrete  can  be  represented  by  a  transfer  function  G(s)  which  is  a  function 
of  the  complex  variable  s»  The  transfer  function  is  unique  for  the  material 

at  a  given  temperature.   It  is  possible  to  obtain  this  function  experimentally 
from  a  series  of  sinusoidal  load  tests  on  the  material. 

2.  The  transfer  function  for  an  asphaltic  concrete  is  of  the  form 

A(s  +  a  )(s  +  a  )(s  +  a  } 

G^S)  =   (s  +  bjfs  +  b.)(s  +  bfyis  +  b-  )  v'12) 

1      2      3  '     **■ 

where  A  is  a  constant,  and  a.  "s  and  b„  's  are  roots  of  the  numerator  and 
denominator,  respectively. 

3=  The  roots  of  the  denominator  of  G(s),  namely  b,,  b  ,  b  and  b,  ,  are 
negative,  real  and  distinct,  thus  indicating  that  bituminous  concrete  behaves 
as  an  overdamped  system. 

k.     Parameters  obtained  in  the  transfer  function  for  bituminous  concrete 
are  believed  to  be  better  indicators  of  material  performance  than  those 
commonly  used,  such  as,  Young's  modulus,  Poisson's  ratio,  etc.  which  change 
with  rate  and  time  of  loading. 

5.  Temperature  is  the  one  single  factor  which  has  the  greatest  effect 
on  the  transfer  function  of  asphaltic  concrete.   Increase  in  temperature 
increases  the  value  of  the  constant  A  in  the  transfer  function  equation, 

6.  The  transfer  function  is  a  powerful  tool  useful  in  predicting  the 
displacement  of  asphaltic  concrete  under  an  applied  load,  dynamic  or  static. 
By  treating  the  static  load  as  a  step  function  of  time,  the  resulting  dis- 
placement can  be.  calculated  by  means  of  the  transfer  function  derived  experi- 
mentally from  the  dynamic  test.  The  excellent  agreement  between  the  calculated 
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and  measured  values  of  the  displacement  in  this  investigation  validates  the 

concept  that  the  transfer  function  represents  a  material  property  which  is 

independent  of  the  type  of  load  input. 

7-   Through  the  use  of  the  transfer  function  and  without  assuming  any 

spring-dashpot  model,  it  is  possible  to  represent  the  time-dependent  behavior 

of  asphaltic  concrete  by  a  fourth  order  linear  differential  equation  with 

constant  coefficients.   The  coefficients  can  be  computed  from  the  roots  of 

the  denominator  and  numerator  of  the  transfer  function. 
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